Detection of a Companion Lens Galaxy using the Mid- infrared 
Flux Ratios of the Gravitationally Lensed Quasar H1413+117 
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ABSTRACT 

We present the first resolved mid-IR (ll^m) observations of the four-image quasar lens 
H1413+117 using the Michelle camera on Gemini North. All previous observations (optical, 
near-IR, and radio) of this lens show a "flux anomaly," where the image flux ratios cannot be 
explained by a simple, central lens galaxy. We attempt to reproduce the mid-IR flux ratios, which 
are insensitive to extinction and microlensing, by modeling the main lens as a singular isothermal 
ellipsoid. This model fails to reproduce the flux ratios. However, we can explain the flux ratios 
simply by adding to the model a nearby galaxy detected in the H-band by HST/NICM0S-NIC2. 
This perturbing galaxy lies 4'.'0 from the main lens and it has a critical radius of 0'.'63 ± 0'.'02 
which is similar to that of the main lens, as expected from their similar H-band fluxes. More 
remarkably, this galaxy is not required to obtain a good fit to the system astrometry, so this 
represents the first clear detection of an object through its effect on the image fluxes of a gravi- 
tational lens. This is a parallel to the detections of visible satellites from astrometric anomalies, 
and provides a proof of the concept of searching for substructure in galaxies using anomalous 
flux ratios. 

Subject headings: galaxies: structure, gravitational lensing 



1. Introduction 

Ever since cosmological simulations of cold dark 
matter (CDM) indicated there should be many 
more satellite galaxies around the Milky Way than 
detected, it has been an importan t goal to try 
to detect these "missing satellites" ( Klypin et al.l 
19991 iMoore et alll999fh In the case of the Milky 
Way, the Sloan Digital Sky Survey has steadily 
found additional, faint satellites, but with nowhere 
near the expected abundances of subhaloe s (e.g. 
Willman et aDl2005t iBelokurov et alJl2007h . The 
simplest solution in the context of CDM models 
is to suppress star formation in low-mass satel- 
lites, probably through heating and baryonic mass 
loss as the universe re- ionize s to leave a popu- 
lation of dark satellites (e.g. Klvpin et al. 1 19991 : 



Bullock et aDl2000h . Locally there is some hope of 
finding these dark satellites through gamma rays 
emitted by dark matter annihilation, but the like- 
lihood of detec tion depends heavily on the CDM 
properties (e.g. Strigari et al.ll2008h . 

Currently the only other technique able to 
search for dark substructures is gravitational 
lensing, where satellites affect both the im- 
age positions a nd fluxes. In s ome lenses (e.g. 



MG0414+0534 iRos et alJ \200d MG2016+112, 
' KochanedliuOa More et all 120091: HE0435-1223, 
Kochanek et al. 20061 ) satellite galaxies of the pri- 
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mary lens can be detecte d through their eff e cts on 
ima ge positions. But , as Mao fe Schneiderl (|l998f) 
and iMetcalf fe Zhao! (2002) emphasized, image 
fluxes are sensitive to perturbations from very 
low-mass satellites, and we observe many lenses 
with "flux ratio anomalies" where the relative 
image brightnesses cannot b e explained by sim 
pie, central lens galaxies (see Evans fc Witti | 200; 
Kochanek fc Dalall2004HCongdon fc KeetonhoO 



Yoo et al J l2006al lhh. Image fluxes can be altered 
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by two physical effects: granularity in the grav- 
itat ional field created by stars ("microlensing," 
see IWambsganssi 120061 ) or satellites ( "millilens- 
ing"), and propagation effects in the ISM of the 
lens galaxy (scattering/extinction). At radio fre- 
quencies, where absorption is negligible and the 
sources are too large to be significantly affected by 
microlensing, a simple explanation for the anoma- 
lous flux ratios in four-image radio lenses is per- 
turbations due to the expected CDM substruc- 
ture. N othing has altered this basic conclusion 
see e.g.lDalal fc Kochane^2002tlMetcalf fc Zha. 
20021 iKoopmans et all 120031: lEvans fc Witt! 12003 
Congdon fc Keetonll2005HYoo et al.ll2006at Dobler 
& Keeton 2006), although there is debate about 
whether som e of the perturbers may be along the 
line of sight dChen et al. 2003; Wambsg anss et al 



20051 iMetcalfl 120051 ) and about the fraction of 



substructures that survive a t the typical impact 
parameter of lensed images (|Zentner et al. 2005; 
Madau et al.ll2008HXu et al.ll2009f ). To date, how- 
ever, no visible satellite or companion galaxy has 
been convincingly detected using flux ratios rather 
than astrometry; therefore, we lack a "proof of 
principle" demonstration of the flux ratio ap- 
proach. 

The need to fully disentangle the physical ef- 
fects is illustrated by the fact that one optica l 
anomaly included by iDalal fc Kochanekl ([20021 ) . 
that of PG111 5+080, was shown by the mid-IR 
observations of Chiba et al. ( 20051 ) to be due to 
microlensing rather than substructure. Like ra- 
dio emission, the mid-infrared is also insensitive 
to extinction and microlensing, so it is an ideal 
bandpass in which to search for substructure in 
the lens. Quasar spectra rise longward of 1/im in 
vF v due to a du st "torus" surrounding the AGN 
( Antonucci 1993), and the mid-IR emission of the 
torus is too large (~0.1 mas) to be affected by stel- 
lar microlensing. In addition, the mid-IR source 
should have a variability time scale long compared 
to the lens time delays, leading to little variabil- 
ity in the mid-IR flux ratios (e.g. 
2007h . 

In this p aper we discuss t he lens system 



Poindexter et al 



H1413+H7 (jMagain et alJ Il988h . for which we 



have obtained mid-IR (11/im ) data with the 
Michelle camera dRochd [20041 on Gemini North. 
H1413+117 is a four-image quasar lens with a 
quasar redshift z s = 2.55. The lens galaxy is 



faint and difficult to detect against the glare from 
the bright quasar Ther e are position estimates 
fromlKneib et all (|l998bh and lChantrv fc Magainl 
(|2007|) based on the same HST/NICMOS-NIC2 H- 
band image, but no lens redshift. T here are many 



ubli s hed mode l s of H1413+117 ( Kavser et al 
199(1 lYun et all [19971 iKneib et al 



1998a; Chae 



& Turnshek 1999) using a broad range of mass 
distributions. The models have never had diffi- 
culty matching the image positions, but all have 
failed to reproduce the flux ratios. 

A difficulty for all these previous studies is that 
the i mage fluxes are likely affected by both ex tinc- 
tion (jTurnshek et al.lll997t iFalco et al.lll999h and 
(chromatic?) microlensing (see Ostensen et al 



1997tlAngonin et al.ll99oHchae et al ll200ll ; Hutse 
mekers 1993). Our 11/im data essen tially elimi- 
nate b oth of these complications. For a lCardelli et al 
(Il989h Rv = 3.1 extinction curve, the absorp- 
tion at 11/im compared to that at H-band is 
only Ru/Rh — 0.05, approaching 0.07 only as 
the lens redshi f t app roaches that of the source. 
Turnshek et alJ (1997) estimated a maximum ex- 
tinction for H1413+117 relative to image C of ap- 
proximately 0.1 < Ah < 0.3, depending on the as- 
sumed dust model and redshift, corresponding to 
negligible corrections (An < 2%) at 11/im given 
our measurement uncertainties. In the quasar 
rest frame, 11/im corresponds to 3/im, a wave- 
length where the emission should be dominate d 



by the large dust torus (Sa nders et al. 1989). 



For microlensing by stars in the lens galaxy (at 
ziens = I), the Einstein radius of the stars is 
1O 16 (M/O.25M ) 1 / 2 cm, so the torus radius is 
roughly 100 times larger than the Einstein radius 
of the typi cal star and microlensin g effects will be 
weak (e.g. iRefsdal fc Stabelllll99"l . 

The par ticular model we advance here was con- 
sidered by IChae fc Turnshek! (|l999h . where it had 
difficulty reproducing the optical flux ratios. We 
improve on their results by presenting models 
based on mid-IR flux ratios with their insensitiv- 
ity to extinction and microlensing. In Section[2]we 
present our mid-IR observations of H1413+117. In 
Section[3]we present a simple model of the lens sys- 
tem that reproduces the observed flux ratios based 
on perturbations from a faint companion galaxy. 
This galaxy is not required to obtain a good fit to 
the system astrometry, so this represents the first 
clear detection of an object (albeit a visible ob- 
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ject) through its effect on image fluxes. This is a 
parallel to the detections of visible satellites from 
astrometric anomalies. In Section [4] we discuss our 
results. 

2. Flux measurements and errors 

We observed H1413+117 with the Michelle 
camera on Gemini North at 11.2 microns (F112B21 
filter) on 5 Aug 2005. The observing time was 1/2 
hour, of which 493.6 seconds were spent on source. 
The data were processed with the standard Gem- 
ini pipeline mireduce. Our analysis starts with the 
coadded chop and nod subtracted image. These 
initial images have 16 vertical stripes due to the 
16 readout channels; we corrected for this by sub- 
tracting the median of each stripe, masking the 
region containing the quasar images. This proce- 
dure also removes any residual sky flux. Figure [1] 
shows the resulting image. 

We fit the four lensed images using a 2-D Gaus- 
sian to model the PSF, fixing the relative im- 
age positions to those measured from HST images 



(jTurnshek et al.lll997t h Since the sky is so much 



brighter than the quasar images, we use the stan- 
dard deviation of the background pixels as an es- 
timate of the errors of the flux of each pixel. The 
model has nine free parameters: the fluxes of the 
four images (4), the position of image A (2), the 
standard deviation of the Gaussian PSF along the 
two principal directions (2), and the rotation an- 
gle of the Gaussian PSF. We found the best-fit 
model using non-linear Levenberg-Marquardt op- 
timization. We used a Markov chain Monte Carlo 
(MCMC) with 10 4 steps to compute the parame- 
ter uncertainties, and we have run multiple chains 
to check the convergence of the best-fit parame- 
ters and their uncertainties. We also estimated 
the errors by simply computing the goodness of 
fit as a function of each image flux, marginalizing 
over all other variables. Figure Q] shows the best- 
fit model and its residuals, and Figure shows 
the MCMC and Ax 2 estimate for the fluxes (in 
arbitrary units) and their uncertainties. We cali- 
brated the flux measurements using observations 
of the standard star 20 Boo and obtained a total 
flux for H1413+117 consistent with the value from 
Aussel et all (|l998l ) of - 30 mJy. The image flux 
ratios (relative to image A) are listed in Table [TJ 
Table [1] also compares our flux ratios to the 



near-IR values derived from the HST/NIC MOS- 
NIC2 H-band im age bvlKneib et al. (1998b, whose 
values agree with lChantry fc Ma gain 2007]) as wel l 
as the values adopted bv lChae fc Turnsheki ([l999) 
based on optical/UV HST photometry with cor- 
rections for extinction (assuming SMC-like dust 
at redshift 2.55) and for microlcnsing in image 
D (based on the li ne-to-continuum flux ratios of 
iHutsemekers! Il993l ) . Our flux ratios are s ignifi- 
cantly inconsistent with Kneib et al. ( 1998bl ). with 
a x 2 of 8.2 for the two m easurements agreeing, an d 
mildly inconsistent with Chae fc Turnsheki ([1999) , 
with ax 2 of 2.1. Th e differences are dom inated by 
images D and B for Kneib et al. I (Il998bh and im- 
ages C and B for lChae fc Turnsheki (|l999h . These 
characterizations of the differences include the un- 
certainties in both measurements, b ut given our 
uncertainties, the central values of iKneib et al 



(|l998bh and IChae fc Turnsheki il999j) are ruled 
out with x 2 = 9.4 and 6.4, respectively. 
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Fig. 2. — Probability distribution from the MCMC 
method for the H1413+117 image fluxes (which are in ar- 
bitrary units). The mean and lcr error bars determined 
from the best-fit Gaussian curve are shown by the vertical 
dashed lines. These are compared to the best-fit flux value 
and 68% confidence limits (or parameter values at which 
Ax 2 = f from the minimum x 2 ) that result from the x 2 
minimization technique (see text), which are shown by ver- 
tical dotted lines. The Gaussian width (<r) is listed in the 
top right corner and is used in our final results. 



3. Models 

With the intrinsic flux ratios in hand, where 
by "intrinsic" we mean independent of microlcns- 
ing and ISM absorption, we can now explore 
the origin of any flux ratio anomalies. We 
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Fig. 1. — The left, middle, and right panels show the 11 micron data, the best-fit model based on 2-D Gaussians, and the 
model residuals, respectively. 



Table 1 

H1413+117 Image Flux Ratios. 



Image mid-IR Kneib et al. Chae & Turnshck 

A = 1.00 = 1.00 = 1.00 

B 0.84 ± 0.07 0.92 ± 0.04 0.96 ± 0.10 

C 0.72 ± 0.07 0.71 ± 0.03 0.59 ± 0.10 

D 0.40 ± 0.06 0.57 ± 0.02 0.40 ± 0.06 



Note. — The left column shows our best-fit mid-IR flux ratios and 
1-sigma MCMC errors. The other c olumns show corres ponding ratios 
from HST/NICMOS-NIC2 H-band l|Kneib et al.lll998bT l and HST op- 
tical /U V measurements afte r corrections for extinction and microlens- 
ing l lChae fc Turnshelj[T99Sh . 



mode l the system with LENSMODEL0 (jKeeton 



200 lh u sing the HST WFPC/ WFPC2 image po- 
sitions (jTurnshek et alj Il997t see Tabled]). We 



considered three estimates of the lens position, 
all derived from the same HST/NICMOS-NIC2 
F160W (H-band) im age: that from iKneib et al.l 
(|l998bf ). that from I Chantry fc Magainl (|2007L 
and our own estimate. The three estimates (see 
Table [5]) differ, although our estimate is more 
consistent with I Chantry &: Magainl (|2007l ). For 
our standard model we use a relatively loose 
constraint on the lens position of (Aa, AS) = 
(-C/14 ± 0'.'02,0'.'56 ± / /02), which encompasses 
these results and the likely systematic problems in 
determining the position of the faint lens. We also 



include weak priors on the ellipticity of the main 
lens (1 — b/a = 0.0 ± 0.5) and the external shear 
(7 = 0.05 ± 0.05). The main lens galaxy (Gl) 
was initially modeled as a singular isothermal el- 
lipsoid combined with an external shear. This is 
the simplest plausible model for a lens. There is 
considerable evidence in favor of the isothermal 
profile, and properties of 4-image systems other 
than time delays are v ery insensitive t o the radial 
mass distribution (see lKochanekll2006l ). 

This simple model has no difficulty fitting the 
image positions, but cannot reproduce the mid- 
IR flux ratios. We find similar results when us- 

1998bh 



ing the lens positions of 



Chantry fc Magainl (|2007l ). If we use the very 



Kneib et al 



or 



1 http: / /redfive. rutgers.edu/~keeton/gravlens/ 



small position errors for Gl from the latter study 
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Table 2 

H1413+117 Image and Lens Positions. 



Component 


Aa (") 


AS (") 


A 


0.000 ±0.003 


0.000 ±0.003 


B 


-0.744 ±0.003 


0.168 ±0.003 


C 


0.492 ±0.003 


0.713 ±0.003 


D 


-0.354 ±0.003 


1.040 ±0.003 


Gla 


-0.112 ±0.02 


0.503 ±0.02 


Gib 


-0.1365 ±0.0024 


0.5887 ±0.0035 


Glc 


-0.142 ±0.02 


0.561 ±0.02 



Note. — Negative RA values are Eastward of Image A. Im- 
age positions are taken from iTurnshe k et al.| lll997f). All 
three lens position estima tes, Gla IIKneib et al.lll998bh , Gib 
l lChantrv fc M again 20071). and Glc (our analysis) are based 
on the same HST/NICMOS-NIC2 F160W image. 



(0'.'004), the fit is dominated by the lens position 
(Xgi — 40), but the fluxes are still poorly fit. We 
feel that our error estimates for the position of 
the lens galaxy (found by bootstrap resampling 
the dithered sub-images and refitting the model) 
are likely to be more realistic. For our standard 
model, we find % 2 = 21.4 for DOF = 6 degrees 
of freedom. The models have 15 constraints (the 
positions of 4 lensed images and the lens galaxy, 
the 3 flux ratios, plus the priors, albeit weak, on 
the ellipticity e and shear 7) and 9 free parameters 
(the position, mass scale bai, ellipticity e, position 
angle 9 e of Gl, the amplitude 7 and position angle 
# 7 of the external shear, plus the source position). 
Table[3]presents the best-fit model parameters and 
a breakdown of the contributions to the goodness 
of fit. We examined whether the problem was spe- 
cific to a particular image by sequentially broad- 
ening the flux error bars for each image. We find 
improvements in the x 2 of 11.5, 0.2, 1.3, and 8.8 
when we relax the flux constraints on images A, B, 
C, and D, respectively, suggesting that the prob- 
lem lies in the fluxes for the saddle point images 
A and D. 

If we now consider a sequence of models adding 
a singular isothermal sphere (SIS) as a second lens 
on a grid of positions, we find the largest improve- 
ment in x 2 near the position of a galaxy (G2) 
lying at (— 1'/87, +4'/15) from image A (see Fig- 



ure El) . This galaxy is object #14 in iKneib et al 
(|l998al ) and has an H-band flux of H= 20.0 ± 0.2 
mag that is similar to their estimate for the flux 
of Gl (20.6 ± 0.5 mag) . This was also noted in 
Chae fc Turnshekl (jl999h . where in their models 



the flux of image D was still poorly fit. If we sim- 
ply add an SIS at the observed position of G2, 
we significantly improve the fit to % 2 = 4.9 for 
DOF = 5. We also find better fits to the flux ra- 
tios wh en using the lens posit i ons of IKneib et al 



(|l998bh or lChantrv & Magainl (|2007h . G2 is close 



enough to the lens (at tq2 — 4'.'0) that it con- 
tributes more to the model than just an external 
shear and convergence. If we attribute the shear 
in our first model to G2, then we know that G2 
induces a shear of order 7^2 = ^G2/(2rc2) — 0.14 
and that it will produce higher order perturba- 
tions beyond the second order shear and conver- 
gence of order 7G2^Gi/ r G2 — 0.02. If we com- 
pare this to the eigenvalues of the magnifications, 
(A_,A + ) = (-0.15, 1.1), (0.20, 0.86), (0.20, 0.92), 
and (—0.24, 1.12) for images A through D respec- 
tively in our first model, we see that these higher 
order terms can significantly alter the fluxes of the 
images. This then leads to the A% 2 = 16.5 im- 
provement in the fit. We also note that the orien- 
tation 6> 7 of the external shear in our first model is 
nearly aligned with the direction of G2. The best- 
fit parameters of the two-galaxy model are also 
reported in Table [3] Note that the critical radii 
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Fig. 3. — The left panel shows the A\ 2 resulting from varying the position of a second lens galaxy. The right panel is a 
similar grid focused on the main lens region. Only limited effort was expended to force the convergence of recalcitrant pixels. 



of Gl and G2, b G1 = 0"66 and b G2 = 0"63, are 
very similar, as we would expect from their simi- 
lar fluxes. If we vary the mass of G2, we find it is 
well-constrained with b G2 — 0'.'63 ± 0'.'02, despite 
the freedom to simultaneously adjust the external 
shear. 

We investigated whether we had any sensitivity 
to the mass distributions of Gl and G2. First, we 
replaced G2 by a point mass, and found no signif- 
icant change in the goodness of fit (A^ 2 = —0.9). 
This is not surprising, as the changes can only 
come from the difference in mass within the annu- 
lus about G2 that covers the lensed images. Sec- 
ond, we replaced Gl by an ellipsoidal pseudo-Jaffe 
model (p oc r~ 2 (r 2 + a 2 ) -1 ) with the halo trun- 
cated at half the Einstein radius of the SIE model 
(a = 0'.'33). This has a series of more signifi- 
cant effects. In order for Gl to produce the same 
quadrupole moment at the Einstein ring, it must 
have higher ellipticity as it becomes more centrally 
concentrated. Moreover, the required shear and 
ellipticity will scale with the surface den sity at 
the ring, roughly as 1 — n (|Kochanekll200i ). This 
means that a centrally concentrated model for Gl 
will require a flatter model for Gl and a more mas- 
sive G2, leading to changes in the balance between 
the second order terms (shear) and the higher or- 
der terms that could be more easily detected than 
when changing the radial structure of G2. While 
we generally observe these effects when we use the 



truncated model with a = 0'.'33, the increase to 
X 2 = 5.7 (Ax 2 = 0.8) indicates that we also lack 
constraints that can distinguish between extended 
and truncated halos for Gl. 

Interestingly, the orientation of the external 
shear in the two-galaxy model, 6 7 — 50? 1, is 
nearly in the d irection of another galaxy in the 
H-band image ( Kneib et al. 1998a). namely ob- 



ject #1. This suggests a more general analysis 
of these neighbors. Figure 0] shows positions of 
all the nearby objects with significant H-band de- 
tections (objects H 2, 1, 2, 3, 8, 11, and 14 in 
Kneib et al.lll998bh . If we assume these are all in 



the same group as Gl and scale the Einstein radii 
of these objects by their H-band fluxes assuming a 
standard Faber- Jackson relation (L oc <r 4 ), we can 
estimate the deflection scale of each galaxy near 
the lens as b n = ba2(L n /La2) 1 ^ 2 , using the fact 
that b oc a 2 for an SIS. Table 2] lists our estimates 
for the critical radii, b n , shear, j n = b n /2r n , and 
the amplitude of any higher order perturbations, 
b n bGi/ r ni f° r ea ch galaxy, where we use our best- 
fit critical radius of 0'.'6 for b G 2- The galaxies in 
Figure U are shown by open circles representing 
their critical radii and filled circles scaled to their 
higher order perturbations, and the orientations 
of the external shear in our single and two-galaxy 
models are indicated with arrows. 

Examining Table [H we see that the next most 
important galaxies are #1 and #8, where adding 
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Table 3 
Modeling Results. 



Model: 


Gl only 


G1+G2 


Gl+G2+#1 


Gl+G2+#l+#8 


y 2 /DOF 

y 2 

Atmage 


21.4/6 


4.9/5 


3.4/4 


2.9/3 


0.3 


0.1 


0.1 


0.1 


^Cflux 


15.3 


2.4 


1.2 


0.7 


Xgi 


2.0 


1.6 


1.8 


1.8 


x 2 ■ 

A,prtor 


3.8 


0.8 


0.3 


0.3 


bci 


0C69 


0'.'66 


/ /63 


Of! 62 


Aa G1 


-0"167 


-0'.'166 


-0'/168 


-0'/168 


AS G1 


0'.'549 


C/556 


0'.'555 


07554 


e 


0.22 


0.26 


0.25 


0.27 


Oe 


-30?2 


-36?5 


-35? 9 


-37? 3 


7 


0.14 


0.087 


0.057 


0.045 


6 j 


38?5 


50? 1 


59? 3 


66? 


bG2 




tf'63 


0'.'56 


0'/65 


Aac2 




= -l'/87 


= -l'/87 


= -l'/87 


A5 G2 




= 4'.' 14 


= 4'.' 14 


= 4'.' 14 


b#i 






0'.'48 


0'.'75 


Aq#i 






= 4'.' 50 


= 4"50 


A5 #1 






= -6'/27 


= -6'/27 










0"036 


Aa #8 








ee 4'/25 


A5 #8 








EE 0'/61 



Note. — Position angles are measured East of North and positions are relative 
to image A (negative RA values are Eastward of Image A). In rows 2-4, the \ 2 is 
broken down into image positions (Xi ma9e )i flux ratios (x}i ux ): Gl position (Xqi), 
and priors on e and 7 (Xp rior )- 



7 



I 1 


i i 1 , 
„ G2=14 

o 


I i i 

o 11 : 


- 


sir 

A 


- 

o 8 - 








LO 3 




o 2 : 


1 N 




O 1 : 


i 


i 


. , . o H2 : 


-5 




Aoc (orcsec) 


5 



Fig. 4. — Locations of the neighbors of H1413+117 de- 
tected in the H-band llKneib et alj|l998af ). Galaxies are 
plotted as open circles with radii equal to to their critical 
radii estimated from H-band magnitudes, and as filled cir- 
cles with sizes scaled to their nonlinear effects listed in the 
last column of Table [4] The double-headed arrows indi- 
cate the direction of the external shear in the single galaxy 
model (67,1) and two-galaxy model ($7,2). 



an object near #1 will also capture the primary 
effects of galaxies #2 and H2. If we add an SIS 
at the position of #1 we find modest improvement 
in the fit, with x 2 = 3.4 for DOF = 4. Adding 
#8 as a third galaxy does not improve the fit be- 
cause its position is roughly orthogonal to the pre- 
ferred shear direction, and we are using a prior of 
7 = 0.05 ± 0.05 on the shear amplitude. When 
adding both #1 and #8, we find x 2 = 2.9 for 
DOF = 3. The results for these 3- and 4-galaxy 
models are summarized in Table [3] Note that the 
discrepancy between the best-fit critical radii in 
our models and the estimated critical radii in Ta- 
ble [4] is not surprising given the uncertainties in 
the Faber-Jackson relation, in the H-band magni- 
tudes, and in our best-fit critical radii. Part of 
the improvement in these multiple-galaxy models 
is due to a better fit to the flux ratios, but much 
is due to the reduced external shear and our shear 
prior. These results suggest we are seeing some 
effects from these other nearby galaxies, but we 
have also reached the limits of what we can mea- 
sure given the available data. 



4. Discussion 

A persistent embarrassment in lens models has 
been the inability to fit image flux ratios. This 
is justified, in many cases correctly, by the com- 
plexities of the physical effects that modify flux 
ratios: absorpt ion by the i nterst ellar medium in 
the lens (e.g. iFalco et all Il999h . intrinsic time 
variability modulated by the time delays, gravi- 
tatio nal microlensing by stars (e.g. Wambsganssl 
2006), and more massive satellites or neighbor- 
ing galaxies. We know from the structure of 
the Einstein ring images of quasar host galax- 
ies or lenses with other additional lensed features 
that the solution does not lie in givin g the pri 



mary le ns an odd angular s t ructure (lYoo et al 



2006 ajlbl: lEvans fc Wittl l2Q03t iKochanek fc Dalai 



20041 : lcongdon fc KeetoiJl2005h . 



The advantage of mid-IR flux ratios is that 
they largely eliminate all of these complications 
other than the perturbing effects of satellites. 
Given typical diffe rential extinctions seen in lenses 
( Falco et al. 1999), the extinction corrections in 
the mid-IR are negligible. Qua sars should show 
little mid-IR variability (e.g 



Poindexter et al 



20071 ). and the emission regio ns are too large t o be 
affected by microlensing fe.g. lAgol et al. 2000h . It 
would, however, be wise to verify these latter two 
points empirically by coarsely monitoring a few 
lensed systems in the mid-IR. 

The remaining physical effect, perturbations 
from satellites (a.k.a. CDM substructure), is thus 
well isolated. To date, the failure of a cen- 
tral lens galaxy model to fit image positions 
has always been traced to a lu minous perturb 



ing object (e. g. MG0414+0534 iRos et al 



MG2016+112, Kochanck 2006; More et al 



HE0435-1223,|KochairekeLal 



2000; 



2009: 



2006). That this is 



true for "astrometric anomalies" is not surprising 
because astrometric perturbations are dominated 
by the high-mass end of the satellite mass func- 
tion which will tend to have formed stars. It is 
also well-established that the flux-ratio anomalies 
can be explai ned by adding low-mass substruc- 
tures (see. e.g.lDal al & Kochanck 2002; Yoo et al 



2006a; Congdon & Keeton 2005; Koopmans et al. 
20031: lEvans fc WittJl2003t iDobler fc Keetonll200d: 
Metcalf fc Zhaol l2002h . and the masses of the 



satellites can be so low that there is no sur- 
prise in a failure to observe the responsible ob- 
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Table 4 
Lens Environment. 



ID 


H n (mag) 


r»(") 


6n(") 


In 




G1(=H1) 


20.6 ±0.5 




0.5 






G2(=14) 


20.0 ±0.2 


4.0 


0.6 


0.08 


0.019 


1 


20.7 ±0.3 


8.3 


0.5 


0.03 


0.003 


2 


22.3 ±0.9 


7.3 


0.2 


0.02 


0.002 


3 


22.0 ±0.5 


8.3 


0.3 


0.02 


0.002 


8 


21.8 ± 0.5 


4.4 


0.3 


0.03 


0.007 


11 


22.3 ± 0.8 


7.3 


0.2 


0.01 


0.002 


H2 


21.6 ± 0.7 


11.0 


0.3 


0.01 


0.001 



Note. — H-band magnitudes (H n ), angular separations from 
Gl, estimated critical radii, contributions to shear, and higher 
order, non-linear perturbatio ns for nearby objec ts detected in 
H-band image of H1413±117 l lKneib et ajj|l998al) . 



ject(s). We should also note that the introduc- 
tion of a visible satellite to explain an astromctric 
anomaly does not necessarily explain any accom- 
panying flux rat io anomaly. For example, while 
Ros et all <l2000h showed that a luminous satellite 
(jSchechter fc Moord[i~993h can explain the image 
positions in MG0414+0534, the model still cannot 
explain the flux ratio of the Aj and A 2 images - a 
problem r ecently confirmed by the mid-IR obser- 



(l2009h . 



vat ions of iMinezaki et al 

Nonetheless, high-mass substructures should 
also produce flux ratio anomalies, and we should 
be able to demonstrate that these can be correctly 
detected and attributed. While our original goal 
was to use the mid-IR flux ratios to search for 
evidence of low-mass, non-luminous substructure 
in the lens HI4I3+117, we instead found that the 
flux ratios could be explained by a nearby galaxy. 
This galaxy, which we refer to as G2, is object 
#14 in the Kneib et al. (1998ab) catalog, and it 
lies T/87 East and 4'.'15 North of image A. It has 
comparable luminosity to the main lens, and it is 
close enough to the lens for the higher order terms 
in its perturbations to the lens potential beyond 
an external shear to significantly modify the im- 
age magnifications. Adding other nearby galaxies 
does not significantly improve the fits, in large 
part because we have a reduced \ 2 near unity af- 
ter adding G2. More importantly, no perturbing 
galaxies are required to explain the astrometry - 
the detection of the neighboring galaxy is purely 



due to the flux ratios. The only similar claim 
is for the four- image lens B2045+265, where the 
anomaly can be explained by turning a relatively 
round visible satellite into a 10:1 flattened mass 
distribut ion running across the most anomalous 
image (jMcKean et all 120071) . While the model 
works to explain the anomaly, the mass distribu- 
tion assigned to the satellite is not very physical. 

Independent of any such debates, mid-IR imag- 
ing provides a simple means of greatly expand- 
ing the numbers of lenses where these issues can 
be studied quantitatively. While there are 10 
known radio lenses with four or more images, 
there are 17 such optical lenses. Moreover, ra- 
dio flux ratios must also be corrected for intrin- 
sic variability modula ted by the time delays (e.g. 



Fassnacht et al. 20021 ). while there is at least an 



expectation that these effects will be minimal in 
the mid-IR. If true, single epoch measurements 
suffice in the mid-IR, while monitoring campaigns 
are required for radio wavelengths. Mid-IR flux 
ratios will also help in models for microlensing, 
since a better mass model can be constructed, and 
the microlensing flux anomalies will be better mea- 
sured. Furthermore, a survey of microlensing com- 
paring the optical to the mid-IR (assuming the op- 
tical can be corrected for extinction) could place 
a constraint on the ratio of dark matter to stel- 
lar matter at the position of t he lensed images 
(|Schechter fc Wambsganssll2002l ). 

Unfortunately, little use has been made of this 
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approach. The re are published measurements 
for Q 2237+0305 (|Agol et al.ll2000UMinezaki et al ' 
2009|), B1422+231 and PG1115+080 flChjbaetal 



20051 ). MG0414+0534 (jMinezaki et alj|2009h . and 
now our measurements for H1413+117. All 
of these measurements are of modest preci- 
sion (5-10%), so there is considerable scope 
both for expansion of the sample and improve- 
ments in the precision. Spitzer/IRAC observa- 
tions (e.g. HE1104-1805. IPoindexter et all [2007 ; 
SDSS1004+4112, Ross et al. in prep), have much 
higher precision but can only resolve the larger 
separation lenses due to the poor resolution of 
Spitzer (FWHM of 1"6 and l'/9 at 3.6 and 8.0/Ltm, 
respectively). Moreover, the highest resolution 
bands will be dominated by more compact ac- 
cretion disk emission in the rest-frame of most 
quasars (3.6/zm would be 1/xm in the rest-frame of 
H1413+117). Thus, further progress will depend 
on ground-based observations until the advent of 
the James Webb Space Telescope (JWST), which 
will provide exquisite mid-IR precision for every 
known gravitational lens. These observations are 
challenging for the current generation of 8m tele- 
scopes, but will become relatively easy with the 
next generation of 30m telescopes. 

This work is based on observations obtained 
at the Gemini Observatory (with Program ID 
GN-2005B-Q-43), which is operated by the As- 
sociation of Universities for Research in Astron- 
omy, Inc., under a cooperative agreement with 
the NSF on behalf of the Gemini partnership: 
the National Science Foundation (United States), 
the Science and Technology Facilities Council 
(United Kingdom), the National Research Coun- 
cil (Canada), CONIC YT (Chile), the Australian 
Research Council (Australia), Ministrio da Cincia 
e Tecnologia (Brazil) and Ministerio de Ciencia, 
Tecnologa e Innovacin Productiva (Argentina). 
This work is based in part on observations made 
with the NASA/ESA Hubble Space Telescope. 
Support for program GO-7495 was provided by 
NASA through a grant from the Space Telescope 
Science Institute, which is operated by AURA, 
Inc., under NASA contract NAS5-2655. EA is par- 
tially supported by National Science Foundation 
CAREER Grant No. 0645416. CSK is supported 
by NSF grant AST-0708082. 
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